High energy spin excitations in BaFe2As2 
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We report neutron scattering measurements of cooperative spin excitations in antiferromagneti- 
cally ordered BaFe2As2, the parent phase of an iron pnictide superconductor. The data extend up 
to ~ 100 meV and show that the spin excitation spectrum is sharp and highly dispersive. By fitting 
the spectrum to a linear spin- wave model we estimate the magnon bandwidth to be in the region of 
0.17 eV. The large characteristic spin fluctuation energy suggests that magnetism could play a role 
in the formation of the superconducting state. 
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One of the greatest challenges presented by the re- 
cently discovered h'on pnictide superconductors^ is to 
identify the electron pairing interaction which permits 
the formation of a superconducting condensate. In con- 
ventional superconductors this interaction is provided 
by the exchange of a phonon. For the iron pnictides, 
however, theoretical calculations^i^ indicate that the 
electron-phonon coupling is too weak to account for the 
observed high critical temperatures. Attention has there- 
fore turned to other types of bosonic excitations which 
could mediate the pairing interaction. 

One such candidate is spin fluctuations'*. In common 
with the layered cuprates, superconductivity in the pnic- 
tides is found in close proximity to parent phases which 
exhibit long-range antiferromagnetic orderi^. However, 
unlike the cuprates, whose magnetic properties are gov- 
erned by strong superexchange interactions between lo- 
calized spin-i moments in a single Cu 'idx2_y2 orbital, 
magnetism in the pnictides is more itinerant in character 
and derives from multiple d orbitals. It may also involve 
a degree of frustration. In magnetically ordered materi- 
als the dominant magnetic excitations are coherent spin 
waves. Wavevector-resolved measurements of the spin- 
wave spectrum by inelastic neutron scattering provide in- 
formation on the fundamental magnetic interactions and 
can also reveal effects due to itinerancy and frustration. 
Such studies on the magnetically ordered parent phases 
of unconventional superconductors like the cuprates and 
iron pnictides are important to establish the characteris- 
tic energy scales of the spin fluctuations and also to pro- 
vide a reference against which changes associated with 
superconductivity can be identified. 

Here we present neutron scattering data on the collec- 
tive spin excitations in antiferromagnetic BaFe2As2. We 
find that the spin excitation spectrum has a very steep 
dispersion within the FeAs layers with a bandwidth in 
the region of 0.17 eV, not much less than that in the 
cuprates. Such a high characteristic energy suggests that 
spin fluctuations are a serious candidate to mediate high 
temperature superconductivity in the iron pnictides. 

The parent phase BaFe2As2 becomes superconduct- 




FIG. 1: (Color onhne) Crystal and magnetic structure of 
BaFe2As2. On the left is the crystal structure with the 
conventional unit cell for the low temperature orthorhombic 
Fmmm structure. The three-dimensional antiferromagnetic 
(AFM) ordering of Fe spins is indicated. On the right is a sin- 
gle layer of Fe spins showing the in-plane AFM order and the 
nearest- and next-nearest-neighbour exchange interactions. 



ing on doping with holes'" or on application of pressure^. 
At Tg = 140 K, BaFe2As2 undergoes a structural transi- 
tion from tetragonal to orthorhombic and simultaneously 
develops three-dimensional long-range antiferromagnetic 
orderiSiiSiiii. On cooling through Tg, the space group 
changes from li/mmm (lattice parameters a — 3.96 A, 
c = 13.0 A) to Fmmm (lattice parameters a = 5.61 A, 
b = 5.57 A, c ~ 12.9 A). The magnetic structure of 
BaFe2As2, shown in Fig.[l] is a coUinear antiferromagnet 
with propagation vector Qaf = (1,0,1). The ordered 
moments on the Fe atoms are of approximate magnitude 
0.9 /iB and point along the orthorhombic a axis. 

Polycrystalline BaFe2As2 was prepared by reacting 
stoichiometric amounts of the elements in a tantalum 
ampoule sealed under argon (800°C for 2 days, then 
900°C for 2 days after regrinding). Phase purity was 
confirmed using X-ray powder diffraction. The neutron 
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scattering experiments were performed on the MERLIN 
chopper spectrometer at the ISIS FacihtyiS,. Approxi- 
mately 8g of the BaFe2As2 powder was sealed inside 
a cylindrical aluminium can mounted in a top-loading 
closed-cycle refrigerator. Spectra were recorded at a tem- 
perature of 7K with four different neutron incident en- 
ergies: E[ = 25meV, 50meV, 200 meV and 400 meV. 
The scattering from a standard vanadium sample was 
used to normalize the spectra and to place them on an 
absolute intensity scale with units mbsr~^ meV~^ f.u.~^, 
where f.u. stands for 'formula unit' (of BaFe2As2). The 
spectra were azimuthally-averaged and transformed onto 
a (Q, energy) grid, where Q = |Q| is the magnitude of 
the neutron scattering vector. The presented intensity 
is the partial differential cross-section d^a/dfldEi multi- 
plied by the factor ki/k^, where fcj and fcf are the initial 
and final neutron wavevectors and Ei is the final neutron 
energy. 

Figure [2^a) illustrates the general features of the data. 
At low energies there is strong diffuse scattering due to 
the elastic peak and scattering from phonons, the lat- 
ter of which increases with Q. The phonon signal drops 
off sharply above 40 meV, which is the upper limit of 
the vibrational density of states^. Two distinct features 
stand out from the phonon signal. One is a narrow pillar 
of scattering at Q = 1.2 A~^, and the second is a plume 
of intensity centred on Q = 2.6 A~^. The latter extends 
in energy to at least 90meV where it disappears out of 
the accessible region of (Q, energy) space. The 1.2 A-i 
feature is followed to lower energies in Fig. [2fb), which 
was obtained with a higher resolution configuration. 

The origin of the 1.2 and 2.6 A~^ features is the 
cooperative spin wave excitations (magnons) associated 
with the antiferromagnetic (AFM) zone centres (1,0, I) 
and (1,2, I). Since the structure is layered we expect only 
a weak variation in the inelastic scattering with I. The 
effect of powder-averaging and resolution-folding makes 
the 2D magnon scattering appear at slightly higher Q 
than the 2D AFM wavevectors (5(1,0) — 1-1 A~^, Q[i,2) — 
2.5 A~^, as observed in Fig.[21 

Figure [3] shows examples of a series of cuts taken 
through the data at different energies. At the higher 
energies the signal is seen to broaden (right panel). This 
is due to dispersion of the spin waves. Below ~15meV 
the magnetic signal decreases in intensity. We fitted the 
cuts taken through Q = 1.2 A^^ with a Gaussian line 
shape on a quadratic background and plot the integrated 
intensities of the fitted peaks at each energy in the in- 
sert to Fig. [3l The data show that the magnetic excita- 
tions are gapped, with no detectable signal below 5 meV. 
However, the gap is not sharp since a sharp gap would 
produce the resolution-broadened step in energy shown 
in Fig. [3l The broadening of the step could be due to 
dispersion of the gap in the c direction and/or the ex- 
istence of two or more gaps in the 5-15 meV range. In 
orthorhombic symmetry two gaps are expected at Qaf 
as a result of magnetic anisotropy which splits the spin 
waves into two non-degenerate branches with predomi- 
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FIG. 2: (Color online) Neutron scattering spectra of 
BaFe2As2. The data were recorded at a temperature of 
7K with incident neutron energies of (a) 200 meV and (b) 
50meV. The pillars of scattering centred near 
and Q = 2.6 A~^ are caused by highly dispersive collective 
excitations associated with the antiferromagnetically ordered 
Fe spins (see Fig. [1} • (b) is a higher resolution image showing 
the magnetic signal emerging from Q = 1.2 A^^. 



nantly in-plane and out-of-plane character, respectively. 

We now compare our data with a linear spin wave 
model for BaFe2As2 based on an effective Heisenberg 
spin Hamiltonian. The suitability of such a model is 
perhaps questionable in view of the itinerant character 
of the magnetisroi^, but at the very least it will provide 
an estimate of the scale of the magnetic interactions. We 
calculated the spectrum by the same method as Yao and 
Carlson^^ but extended the Hamiltonian to include terms 
that represent the single-ion anisotropy: 

H = Y^ J,kS, ■ Sk + J2{K,iS^,),+Ka,iS^y - Si),}. (1) 
O'fc) j 

The first summation is over nearest-neighbour and next- 
nearest-neighbour pairs with each pair counted only once. 
The Jjk are exchange parameters as defined in Fig. (TJ 
and Kab and Kc are in-plane and out-of-plane anisotropy 
constants, respectively. Diagonalisation of equation ([T]) 
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FIG. 3: (Color online) The left and right panels show a se- 
ries of constant-energy cuts through the magnon signals at 
Q = 1.2 and Q = 2.6 A~^, respectively. The data were 
averaged over the energy ranges indicated. Successive cuts 
are displaced vertically for clarity. The symbols represent dif- 
ferent neutron incident energies: 25meV (red circles), 50meV 
(blue squares) , 200 meV (green triangles) , 400 meV (grey in- 
verted triangles). In the left panel the lines are fits to Gaus- 
sian peaks on a sloping background. In the right panel the 
lines are constant-energy cuts (averaged over the same en- 
ergy ranges as the data) through the powder-averaged spin 
wave spectrum calculated with parameters SJ2 ~ 2SJia ~ 
2SJu = 35meV, Jc = 0, SKab = SK^ = 0.042 meV and 
S'eff = 0.28 (see Fig. |4|. The insert shows the integrated in- 
tensities of the magnon scattering measured at Q = 1.2 
as a function of energy. The line indicates the expected step 
in intensity for a clean gap folded with the experimental res- 
olution. 
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Aq = 25'{ Jif,[cos(-^) - 1] + Jia + 2J2 + Jc} 



C = S{Kab~K, 



= 25{Ji,cos(^) + 2J2Cos(^)cos(^) 

-KJcCos(Q • c)}. (3) 
The neutron scattering cross section may be written^^ 
(Pa k{ /7ro\2 
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x5](<5a/3-QaQ/3)5"'^(Q,C^), (4) 

a/3 

where = 72.8 mb, g is the g-factor of iron, f{Q) 

the form factor of iron, exp(— 2VF) is the Debye-Waller 
factor which is close to unity at low temperatures, Qa 
is the a component of a unit vector in the direction of 
Q, and ^"^(Qja;) is the response function describing af} 
spin correlations. Only the transverse correlations {yy 
and zz for BaFe2As2) contribute to the linear spin wave 
cross section and the response functions (per BaFe2As2 
formula unit) for magnon creation are given by 

SyyiQ,Lu) = S,s^^—^j^^{niLu) + l}d[uj-coim 
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{n{uj) + l}6[u - uj2{Q)] 
(5) 



where S'ofr is the effective spin and n(uj) is the boson 
occupation number. In linear spin wave theory SeS — S, 
but we keep them distinct here because in the analysis 
they are obtained essentially independently (see below). 

Because our data do not extend over the full spin wave 
dispersion we were unable to determine Jia, Ju and J2 
independently. However, several author si^iiSiiiSiSS have 
made predictions of effective Heisenberg exchange pa- 
rameters from first-principles electronic structure calcu- 
lations. In some casesi^ii^ii^ J2 is predicted to exceed 
Jia and Jib by about a factor of two, with Jia and Jib 
either both ferromagnetic or both AFM (providing Jia 
and Jib are the same sign the spectrum at low energies is 
not very sensitive to which sign it isi^). Alternativeljs^S, 
Jia 3> Jib and Jia — 2J2. Guided by these predictions 
we fixed the ratios of the exchange parameters to be ei- 
ther (i) Jia = Jib = J2/2, or (ii) Jia = 2J2 = -hJib- 
Since we cannot resolve more than one gap in the data 
we set Kab = Kc so that the in-plane and out-of-plane 
gaps are the same, and we neglected the c-axis coupling, 
which is expected to be much smaller than the in-plane 
coupling. The g-factor was set to 2. 

We computed the powder-averaged spin wave spec- 
trum convoluted with the instrumental resolution and 
fitted it to the experimental data near Q = 2.6 A-i al- 
lowing only SJ2 and Soff to vary. Ses is essentially de- 
termined by the absolute intensity, and SJ2 by the dis- 
persion. Data near Q = 1.2 were excluded as the 
peak widths are dominated by instrument resolution and 
consequently are insensitive to SJ2 ■ We obtained equally 
good fits to the data with both sets of parameter ratios. 
The best-fit parameters are, for case (i) SJ2 = 35±3meV 
and Scs = 0.28 ±0.04, and for case (ii) SJ2 = 18±lmeV 
and S'ofr = 0.54 ± 0.05. The highest resolution data re- 
veals a gap of 7.7 ± 0.2 meV — Fig. [3] (insert) — which 
yields SK^ = SKab = 0.042 ± 0.004 meV [case (i)] and 
0.053±0.005meV [case (ii)]. For case (i) we tried different 
values of J1/J2, but the best-fit value of 5/2 was rela- 
tively insensitive, varying from 33meV (J1/J2 = 0.25) 
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FIG. 4: (Color online) Simulation of the powder-averaged spin 
wave spectrum of BaFe2As2. The simulation covers the same 
Q and energy range as the data in Fig. [2ja) . The parameters 
used for the simulation are given in the caption to Fig. O 



to 46meV (J1/J2 = 1). Aside from the uncertainty in 
J1/J2, the major error in SJ2 is the statistical error on 
the fit. The error in ScB comes from estimates of the 
background and instrumental resolution. 

To illustrate the level of agreement, Fig. |4] shows the 
powder-averaged spin wave spectrum calculated with the 
best-fit parameters for case (i) over the same range of Q 
and energy covered by the measurements in Fig. [U^a). 
The maximum spin wave energy is ~175 meV for both 
parameter sets, well out of range of the present exper- 
iment. We note that the order of magnitude of the 
experimentally-determined exchange parameters is con- 
sistent with theoretical prediction o^^'^^i^^'^" providing S 
is close to unity. The values of gScs [0.56 for case (i), 1.08 



for case (ii)] obtained from our analysis are comparable 
with the measured ordered moment in /xb of 0.8-0.9^^'—. 

Our results show that spin fluctuations in the parent 
phase of the pnictides exist over a wide energy range 
extending up to ^175 meV, not much less than that 
found in the cuprate high temperature superconductors. 
Recent neutron scattering studies on single crystals of 
SrFe2As2^^ and CaFe2As2^'^, though restricted to ener- 
gies below 25meV, also find steep magnon dispersion 
relations. For high Tc it is natural to look for a pair- 
ing boson with a large characteristic energy. The data 
here show that spin fluctuations satisfy this requirement. 
However, if magnetism is to play a role in the supercon- 
ducting state then there must also be a coupling between 
spin fluctuations and the electrons involved in pairing. 
One piece of evidence for this is the recent observation of 
a resonant spin excitation in the superconducting state of 
Bao.6Ko.4Fe2As2^. More generally, the itinerant charac- 
ter of magnetism in the pnictides has been inferred from 
the nesting of electron and hole Fermi surface pockets 
with nesting vector Qap^>^'^>^ and by the observed ef- 
fect of AFM order on the Fermi surface as revealed for 
example by optical spectroscopj«^i^, angle-resolved pho- 
toemission spectroscopy^ and quantum oscillations^^. 

On the other hand, the dynamic magnetic response 
measured here does not show any obvious fingerprints 
of itinerant magnetism, such as damping due to a 
Stoner continuum. It will be interesting to follow the 
spin excitation spectrum to still higher energies where 
itinerant effects generally have the largest influence. 

We thank 1.1. Mazin for helpful comments. This work 
was supported by the Engineering and Physical Sciences 
Research Council of Great Britain. 



* Electronic address: |a.boothroyd@physics.ox.ac.uk 

^ Y. Kamihara et al, J. Am. Chem. Soc. 130, 3296 (2008). 
^ L. Boeri, O. V. Dolgov, and A. A. Golubov, Phys. Rev. 

Lett. 101, 026403 (2008). 
^ K. Haule, J. H. Shim, and G. Kotliar, Phys. Rev. Lett. 

100, 226402 (2008). 

* D. J. Singh and M.-H. Du, Phys. Rev. Lett. 100, 237003 
(2008). 

^ J. Dong et al, Europhys. Lett. 83, 27006 (2008). 
" C. Cruz et al, Nature (London) 453, 899 (2008). 
^ M. Rotter, M. Tegel, and D. Johrendt, Phys. Rev. Lett. 

101, 107006 (2008). 

® P. L. Alireza et aZ.. larXiv:0807n8 96! (unpublished). 

^ M. Rotter et al., Phys. Rev. B 78 020503(R) (2008). 

" Q. Huang et al, arXiv:0806.2776 (unpublished). 

^ Y. Su et ai, arXiv:0807.1743 (unpubhshed) . 

^ R. I. Bewley et al, Physica B 385-386, 1029 (2006). 

^ R. Mittal et al, Phys. Rev. B 78 104514 (2008). 

-'^ I. L Mazin et ai, arXiv:0806.1869 (unpublished). 

^ D. X. Yao and E. W. Carlson, Phys. Rev. B 78 052507 



(2008). 

G. L. Squires, Introduction to the theory of thermal neutron 

scattering (Dover Publications, Inc., 1996). 

T. Yildirim, Phys. Rev. Lett. 101, 057010 (2008). 

Q. Si and E. Abrahams, Phys. Rev. Lett. 101, 076401 

(2008). 

F. Ma, Z-Y. Lu, and T. Xiang, arXiv:0804.3370l (unpub- 
lished). 

Z. P. Yin et ai, Phys. Rev. Lett. 101, 047001 (2008). 
J. Zhao et al, arXiv:0808.2455 (unpublished). 
R. J. McQueeney et al, arXiv:0809.1410 (unpublished). 
A. D. Christianson et al, arXiv:0807. 3932.. (unpublished). 
V. Cvetkovic and Z. Tesanovic, larXiv:0804.4678l (unpub- 
lished) . 

F. Ma and Z-Y. Lu, Phys. Rev. B 78, 033111 (2008). 
K. Kurokiei ai, Phys. Rev. Lett. 101, 087004 (2008). 
W. Z. Hu, et al, arXiv:0806.2652 (unpublished). 
L. X. Yang et ai, arXiv: 0806. 2627 (unpublished). 
S. E. Sebastian et ai. arXiv:0806.472 6 (unpublished). 



